Two photon ionization of condensate atoms 
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The efficient photoionization of a Bose-Einstein condensate requires the creation of ions with the 
smallest possible transfer of atoms from the condensate into the thermal phase. The spontaneous 
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decay from excited states into atomic states with momentum different from the initial one should 
be reduced. We investigate theoretically the two-photon ionization of a rubidium condensate using 
near resonant excitation to the 6P state and second photon at 421 nm or 1002 nm into the contin- 
uum. Different ionization schemes with coherent control of the first excitation and reduction of the 
spontaneous decay are presented. 
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PACS numbers: 03.75.Nt,32.80.Rm,42.50.Hz 
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^ ' I. INTRODUCTION 

The realization of Bose-Einstein condensates (BEC) of alkali atom vapors has attracted much interest into new 
aspects of photon-matter interaction arising from the quantum nature of the atomic sample. Recently attention has 
been paid to the photoionization of a BEC by monochromatic laser light Q, 0, 0, 0, 0] ■ The products of the pho- 
toionization process (electrons and ions) obey to Fermi-Dirac statistics. Due to the quantum nature of the initial 
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atomic target and the narrow spectral width of the laser ionization source, the occupation number of the electron/ion 
. i-h . final states may become close to unity, especially for a laser excitation close to threshold. If this regime is reached 
the ionization rate could be slowed down by the Pauli blockade, with the occupation of ionized states determined by 
the balance between the laser ionization of the condensate and the rate of escape of the ionization products from the 
condensed system. Moreover, the electrostatic interaction between the ions and neutral atoms in the condensate is 
expected to originate localized deformations in the condensate density distribution. Either an increase or a decrease 
in the neutral atom spatial density, depending on the sign of the atom- ion phase shift, could be originated in small 
regions surrounding the ions @ . 

The experimental study of these processes requires the production of a large ionic concentration within a conden- 
\f} ' sate on time shorter than the timescale on which the ions escape from the system, of the order of hundred ns for 
typical condensate samples, without simultaneously destroying the condensate. Instead, in the two-photon ionization 
experiments reported in refs. 0- [J the photoionization laser produced a significant transfer of atoms from the 
condensate into the thermal cloud. Those experimental investigations demonstrated that the rubidium two-photon 
ionization was not performed in efficient way, i.e. without producing a large depletion of the condensate atoms with 
their transformation into thermal atoms. Indeed, in typical two-photon ionization experiments the ionization is in- 
. creased only at the expenses of a large production of thermal atoms. 
. ' The present work reports a theoretical investigation of the efficiency realized in the ionization of atomic rubidium. 

This work addresses a target which is different from those associated with the typical investigations of multiphoton 
^ ^h * processes. Our aim is not to increase the multiphoton ionization probability per se. In contrast, the production of 
ions is compared to optical pumping rate, produced by the spontaneous emission, transferring condensate atoms into 
the thermal cloud. The ionization efficiency will be determined by the decrease in the loss towards atomic states 
initially not occupied, i.e. a decrease of optical pumping into sink thermal states, and by the increase in the number 
of ions. 

Our condensate ionization scheme takes place in a three-level cascade scheme, the final one being the atomic contin- 
uum. The condensate losses are produced by spontaneous emission from the intermediate excited state. An efficient 
ionization is realized for an atomic transfer into the continuum with a reduced real excitation of the intermediate 
level, whence a reduced role of spontaneous emission processes from that intermediate level. The task of reducing the 
transfer towards sink levels is usually reached applying the STIRAP technique for the coherent transfer between the 
initial and final levels of a three level scheme, as analyzed in [ZLIg, and for transfers to and from a continuum as in 
111 [HHm An efficient STIRAP transfer, which relies on the coherent coupling of the initial and the final state, can 
be realized only using very short laser pulses, i. e. shorter than the characteristic decoherence time of the system. 



•present address: National Institute of Standards and Technology, Gaithersburg, MD 20899, USA. 



2 



In the case of photoionization, the use of laser pulses shorter than the typical electron-ion decoherence times would 
prevent from creating the narrow kinetic energy distribution for the electrons and ions required for the Pauli blockade. 
Instead we will analyze here the two-photon ionization produced by narrow band cw radiation. 

For long interaction times the photoionization and optical pumping processes can be described through a simple 
rate equation approach. At shorter interaction times, comparable to or smaller than the spontaneous emission life- 
time from the first excited state, the determination of the photoionization efficiency requires the solution of density 
matrix equations. Our solution of the densi ty m atrix equations searching for a regime of efficient ionization has strong 
connections with the investigations of refs. |12U13| for the lineshapes of atomic excitation and ionization in a cascade 
three-level system under two-color excitation as a function of the time delay in the application of the two lasers. 

Section II describes the characteristics of our three-level scheme. Section III analyzes the case of simultaneous 
application of the excitation and photoionization lasers for long interaction times on the basis of the rate equation 
solution. In Section IV we analyze the ionization efficiency for different temporal shapes of the excitation and ioniza- 
tion laser pulses using the density matrix approach. The lengths of the laser pulses and their time delay represent the 
experimental parameters controlling the ionization efficiency. Section V analyzes ionization based on a rapid adiabatic 
transfer. 




FIG. 1: Three-level scheme for the excitation from the condensate ground state to the condensate ionic state through an 
intermediate |e) state. The thermal sink state \s) is occupied by the spontaneous emission decay from the |e) state. 



II. ATOMIC LEVEL AND HAMILTONIAN 

We investigate the rubidium two-photon ionization scheme from the 5Si/2 ground state to the continuum with 
intermediate resonant state the 6P 3 / 2 state at 2.94 eV above the 5S!/ 2 ground state as in the experiment of ref. |4|. 
The two-photon ionization was produced either by two photons a 421 nm or by a two color process, one photon at 421 
nm and a second one at 1002 nm. The photoionization cross sections for absorption of one 421 nm or one 1002 nm 
photon from the 6P state have a small difference, not playing an important role for the present investigation of the 
photoionization efficiency. Instead a larger intensity was available at the 1002 nm wavelength, and the contribution 
to ionization of the two color process was much larger than the contribution of the one-color photoionization by two 
421 nm photons. 

The level scheme of Fig. 1 schematizes the two-photon ionization process from the ground state |<?) of the condensate 
to the condensate ionic state \i) with a near resonant excitation of the intermediate level e). The population of the 
intermediate state |e) may decay back into the ground state through emission of spontaneous emission photons. For 
a sample composed by a Bosc-Einstcin condensate with atoms having a momentum p = 0, the spontaneous emission 
produces ground state atoms with an atomic momentum p different from zero. These atoms in the p / ground 
state constitute a thermal component separate from the condensate, to be described through a sink state \s) whose 
atomic properties are equivalent to those of the \g) state. The |s) state is filled by the spontaneous emission with rate 
r from the |e) state. The interaction of the \s) state with the photoionization radiation is the same one of the |<?) 
state. However the photoionization from the \s) state produces ions having an initial momentum p =/= and weakly 
interacting with the condensate. Therefore that photoionization process is not relevant for our target of electrons/ions 
with a small kinetic energyQ and it will not be considered into our analysis. 

A laser with angular frequency ui, Rabi frequency f^i and detuning 5 = u)\ — Lu cg acts on the \g) — > |e) transition, 
while a laser with angular frequency 0J2 and intensity I2 acts on the |e) — > \i) transition. The intensity of the first 
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FIG. 2: Comparison between the pertubative solution of Eqs.JHJl (solid lines) and the numerical solution of Eqs. 111-")! (dashed 
lines) for the laser pulse shape shown as dotted line in (a). In (a) fraction p ca , in (b) sum of ionized fraction pa and of sink-state 
fraction p EE . Laser parameters Ji = 0.66 Wcm~ 2 , corresponding to Qi/2tt = 5.0 MHz, I2 = 2.0 x 10 3 Wcm -2 . 
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laser is I\ — 24"^ The ionization from the excited state |e) is described by the Ui ionization cross section for the 
laser with intensity li. For rubidium the atomic parameters are T = 8.9 x 10 6 s -1 , o\ = 4.7 x 10~ 24 m~ 2 for the 421 
nm radiation and 02 = 16.4 x 10~ 24 m~ 2 for the 1002 nm radiation|5|. 

The evolution of the atomic system interacting with the laser radiation is given by the Optical Bloch Equations for 
the atomic density matrix [bH UJl on the Hilbert space spanned by the four atomic states introduced above 

(1) 

- (aih + (T2I2) pee, (2) 

(3) 

-y (Pee ~ Pgg), (4) 

(5) 

with the additional condition pgg + p cc + Pa + Pss = 1, where p eg denotes the interaction representation of p eg . We 
have not included into the density matrix equations the excitation and ionization of the sink state |s) by the two 
lasers, because those processes do not modify the number of atoms lost or ionized from the condensate phase. The 
numerical results reported in the following are valid not only for the simple three- level scheme of Fig. 1, but also for 
the real rubidium level scheme, with the excited state occupation spontaneously decaying to other intermediate states 
before reaching the final sink state[f|- Both the ions and the losses into the sink states deriving from the condensate 
component are equivalent in the two cases. 

The efficiency of the ionization will be measured through the ratio r between the atomic population in the sink 
state and the atomic population in the ionic state, after the end of the laser interaction 

r= /o > ( ^, (6) 
Jo Pn(t)dt 

where the time f is longer than the laser pulse length r in order to include also the spontaneous emission losses 
towards the sink state at the end of the laser-atom interaction. Notice that in this analysis a large (small) value of r 
corresponds to a low (high) ionization efficiency. 

For any temporal shape of the applied laser pulses, a key parameter for the photoionization efficiency is the laser 
detuning S from the \g) — > |e) transition. For decreasing detunings the two- photoionization process becomes more 
efficient owing to the large excitation of the intermediate |e) state; at the same time, however, also the spontaneous 
emission becomes more efficient. As a consequence, for decreasing detunings the loss rate towards the \s) state 
usually increases faster than the ionization rate, resulting in the decrease of efficiency for the condensate atoms 
photoionization. 



4 



III. LONG LASER PULSES 



During the simultaneous application of the two lasers, the atomic population decays from the excited state |e) and 
is lost from the system with total rate equal to IY = T-\-aiI\ + 02/2- Therefore, at every time t the total population in 
the \g) and in the |e) states decays with a time dependent loss rate equal to r L( o ee (i). In the case in which T, o-j/j <C fii 
the ionized fraction of population and the fraction lost by spontaneous emission can be estimated by treating the 
losses as a perturbation of the steady state solution for the Eas. (|ll5|) restricted to the [|<?}, |e)] two-level system. This 
perturbative solution is 
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Supposing the two lasers applied for a time r the ionic and sink state occupations may be written as 

<j\I\ + CT2/2 
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Therefore all the population is lost from the three-level system if the duration t of the pulses is much longer than tl. 
In this limiting case, with no atom remaining in the \g) state, the efficiency of the process is given by 



Km = /. 



(10) 



with the / parameter defined by 



o\I\ + 0-2/2 



For instance, at laser intensity I\ = 0.66W cm 2 corresponding to a Rabi frequency fli/(2ir) 
intensity I2 = 2000 W cm -2 , the / parameter is equal to 53.8. 



(11) 
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IV. PULSED EXCITATION 

As demonstrated in refs [T2I IT^ | , under the interaction with coherent pulsed radiation the atomic excitation presents 
different properties than in measurements performed under the cw excitation regime. In the regime of pulsed exci- 
tation, the time evolution of populations and coherences becomes very important. For calculating the ionization 
efficiency it is necessary to perform the numerical solution of the density matrix equations lllol) . In fact, it should 
be taken into account that the a ctual oc cupation p cc of the excited state is characterized by Rabi oscillations with 
effective Rabi frequency ttf* — \/fl 2 + S 2 . Because the ionization depends on the occupation of that state, the atomic 
ionized fraction increases in time with oscillations at the same frequency. 

Fig. El compares the exact solution of the density matrix equations and the perturbation solution given by Eqs.ljHJ). 
The average ionized and the lost fraction are well described by the perturbative solution. However, due to the oscilla- 
tory evolution of the excited state population, the numerical solution of Eqs. Illl">l) is required in order to obtain the 
correct value of the ionization efficiency parameter r. 

Whenever r < tl, while no more ionization occurs after the duration r of the pulse and the total ionized fraction 
remains equal to Pii(r), a finite fraction of population p cc (t) remains in the excited state at the end of the atom- laser 
interaction and decays by spontaneous emission into the sink state. Therefore for a pulse of duration r the ionization 
efficiency ratio of Eq. is given by 

. . _ Jq Tp ee (t)dt + p cc (t) 
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FIG. 3: In (a) ionized fraction pa and in (b), continuous line, parameter r of the efficiency ratio versus the blue laser detuning 5 
calculated through the numerical solution of the density matrix equations. The dashed line reports the limiting value given by 
the / parameter of Eq. 1111 Laser parameters I\ — 0.66 Wcm~ 2 , corresponding to Qi/2tt = 5.0 MHz, I2 = 2.0 x 10 3 Wcm -2 , 
interaction time r = 0.2/ ps. 



In general r > / and the maximum efficiency of the ionization process corresponds to the minimum value = f. 

The oscillating behavior of the ionization process appears also measuring the ionization as function of the detuning 
S. Such a behavior is shown by the results of Fig. |3|where the ionized fraction pa and the efficiency ratio r are plotted 
as a function of the laser detuning 5 from the |e) excited state, at fixed interaction time r. Owing to the choice 
r = 27r/rZi, Pcc(t) = 0, so that r has a minimum (and the efficiency has a maximum) at zero detuning. At larger 
values of the detunings pa decreases with oscillations, and p ee (T) assumes values dependent on the effective frequency 
of the Rabi oscillations. The role of the Rabi oscillations is more clear on the plot of Fig. ]^b) for the efficiency ratio 
r. 

For fixed pulse duration the oscillating time dependence of the excited state population with detuning-dependent 
effective Rabi frequency il\ s produces an oscillating dependence of the ionized and of the lost fractions, and thus also 
of the ionization efficiency, as a function of the detuning S/2tt. In particular, the large oscillation of the ionization 
efficiency is determined by the oscillating amount of population remaining in the excited state |e) at the end of the 
laser pulse and decaying totally by spontaneous emission into the sink state \s). As a consequece, r assumes oscillating 
values which are always larger than the / value of Eq. Ijllfl . and can be equal to that value only for isolated detunings 
where the effective Rabi frequency fi® ff produces a 2ir pulse for the applied duration of the atom-laser interaction. 

The occupation of the excited state at the end of the laser pulse can be annulled in the regime defined as coherent 
population return^,^]. For the coherent excitation of two- level atomic systems subjected to light pulses of duration 
r and detuning S from resonance, the dynamics of the atomic population is adiabatic if the detuning is significantly 
larger than the Fourier width of the pulse, i.e., if S 3> 1/t. In this case the dressed state occupied at the end of the 
excitation is the same one as before the excitation. For a smooth pulse at constant detuning, that state coincides 
with the atomic ground state, so that even if the excited state is populated during the application of the pulse, all 
the population is adiabatically brought back to the ground state at the end of the process. Therefore no population 
fraction remains in the excited state, p ee (T) = 0. This regime of coherent population return does not produces any loss 
into the sink state at the laser pulse end. Results showing the advantages of this regime in the two-photon ionization 
are reported in Fig. 0] 

We considered several durations of pulses all of them with the shape shown in Fig. Ela). For each duration of the IR 
square pulse, the blue intensity had a Gaussian profile such that the intensity decreased to 1/e 3 of the peak value in 
correspondence of the IR pulse start and end. In Fig. 01(b) and (c) the ionized fraction pa and the parameter r of the 
the ratio between the atomic fraction in the sink state and the ionized fraction are shown, respectively, as a function 
of the laser detuning 5, for several durations of the laser pulse and at a given blue Rabi frequency. In all cases, the 
effects of adiabatic excitation are visible, because for each pulse duration there is a minimum detuning necessary to 
obtain the optimal ratio r. The value of the minimum detuning is only slightly sensitive to the Rabi frequency of the 
transition, because, as pointed out in |l2j, no power broadening occurs in a scheme of coherent population return. The 
limiting value of r reached at large detunings corresponds to the minimum value / calculated at the laser intensities 
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FIG. 4: Ionization with coherent population return. In a) temporal shapes of the blue/IR laser for a 125 ns time length of 
the IR pulse. In b) ionic occupation pa and in c) parameter r describing the efficiency ratio versus laser detuning S/2n, in 
MHz. fii/27r = 5. MHz, h = 1500 W/cm 2 , / = 71.7. Results for different time lengths of the IR pulse, 30, 60, 125 and 250 ns 
respectively, with the blue pulse length scaled as shown in a) and discussed in the text. 

of the numerical analysis. 

V. RAPID ADIABATIC PASSAGE 

Another excitation scheme allowing an efficient and robust transfer of the population from the ground to the excited 
state is the rapid adiabatic passage (RAP) technique 17, 18] . This technique takes advantage of the population 
adiabatic evolution in order to transfer almost the whole population from the ground state to the excited state, 
blocking completely the population return. RAP is produced by a time dependence of the laser detuning, and for an 
efficient transfer the laser detuning changes from a very large and negative value to a very large and positive value, 
or vice versa, during the application of the light. The laser frequency experiences a chirp, and the chirp duration 
should be shorter than the excited state lifetime. Under these circumstances, the lowest energy adiabatic state of the 
system coincides with the ground state before the application of the light and with the excited state at the end of the 
process. When this scheme is applied to the lower atomic transition of the rubidium two-photon ionization, it allows 
to transfer the ground state population to the |6P) state, from where a pulse of IR light, with the time sequence shown 
in Fig. OJleft-hand side), produces the ionization shown in Fig. fright-hand side). Because the IR pulse ionizes the 
excited state atoms and also modifies the energy of the excited state, for any temporal shape of the two laser pulses 
the largest two-photon ionization is obtained for a particular delay of the IR pulse with respect to the blue laser pulse. 
For the pulse shapes of Fig. [5] the maximum ionization was obtained when the center of the IR pulse was delayed by 
0.02 (is with respect to the center of the blue pulse. A linear chirp of the blue detuning was examined, but the chirp 
exact shape had a minor influence on the ionization efficiency, provided that the detuning difference between initial 
and final value was larger than ~100 MHz. Suppressing the oscillatory behavior of the excited population, the RAP 
led to an excited state occupation and an ionized fraction larger than those obtained by pulses of the same shape 
but at constant blue detuning. For the blue laser scanning of Fig. [SJthe ionization fraction was larger by a factor 2 
than produced at constant blue detuning (0.0057 against 0.0027). The RAP technique is thus suitable for producing 
the maximum possible fraction of ionized population within a laser pulse sequence of short duration, during which 
the ratio between the atomic population in the sink state and the atomic population in the ionic state corresponds to 
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FIG. 5: Ionization using RAP from the ground state to the excited state. In left plot the temporal shapes of the blue (dotted 
line) and IR (dashed line) laser pulses and of the blue laser detuning (continuous line) 5/2tt, in MHz, at £li/2ir — 30.7 MHz 
and I2 = 2000 W/cm -2 . In right plot time evolution of the excited state occupation (solid line) and of ionized fraction excited 
state (dashed line). 



the optimal value /. Due to the suppression of the coherent population return, on the other hand, a large fraction of 
population typically remains in the excited state at the end of the laser pulse and decays into the sink state thereafter. 
As a consequence, in RAP at long interaction times the efficiency parameter r, calculated according to Eq. [f)| reaches 
values larger than those associated to the adiabatic regime. For the parameters considered above the efficiency ratio 
r was 165 against the minimum value / = 53.8. 



VI. CONCLUSION 



We have examined a problem to be classified as a quantum control one: reach a specific target, the ionization of 
the condensate, minimizing the losses into a different output channel, the production of thermal atoms. We have 
explored different two-photon ionization schemes, producing different quantum controls for the atomic excited state. 
However the quantum control reached for the two-photon ionization is very poor. In effect the ionization process 
transferring atoms from the excited state to the ionization state is governed by rate equations, whence a classical 
response instead of a quantum one. Our final result is that the maximum ionization efficiency to be reached through 
the quantum control on the first excitation, depends on the probabilities of two processes, the spontaneous decay rate 
and the ionization rate. These rate are not modified within the quantum control of our analysis. 

While the qualitative behavior of the ionized fraction depends only on the properties of the laser pulse exciting 
the first step of the ionization sequence, the absolute value of the ionized fraction depends only on the intensity 
of the IR light coupling the excited state to the continuum. An analogous statement applies also to the efficiency 
ratio between the thermal fraction and the ionized fraction. The criteria we have derived for the optimization of the 
two-photon ionization have a general validity, whereas the absolute numbers we have discussed correspond to the 
intensities available in the analyzed experiment. The fraction of ionized atoms in these conditions is much smaller 
than the number of photons scattered by each atom. The application of the optimization criteria we have derived, 
however, together with IR intensities larger than those we have considered, will allow to reach experimental conditions 
leading to an efficient large photoionization of Bose-Einstein condensates. 
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